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The thermal and nonthermal pion production by sigma 
decay and its relation with chiral symmetry restoration in a 
hot and dense matter are investigated. The nonthermal de- 
cay into pions of sigma mesons which are popularly produced 
in chiral symmetric phase leads to a low-momentum pion en- 
hancement as a possible signature of chiral phase transition 
at finite temperature and density. 
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As is well known, the ultimate goal of relativistic heavy 
ion collisions [1] is to study the deconfinement process in 
moving from a hadron gas to a quark-gluon plasma and 
the chiral transition from the chiral symmetry breaking 
phase to the phase in which it is restored [2]. In recent 
years some signatures of the chiral transition have been 
proposed, such as the excess of pions due to a rapid ther- 
mal and chemical equilibration in the symmetric phase 
[3] , the excess of low energy photon pairs by pion anni- 
hilation in hot and dense medium [4], the dilepton en- 
hancement from leptonic decay of sigma [5], and the 
enhancement of the continuum threshold in the scalar 
channel [6] . Most of the signatures are associated with 
the changes of sigma properties at finite temperatures 
and densities. While in the vacuum there is no obvious 
reason for the presence of sigma mesons in the study of 
chiral properties, there is a definite need [7] for sigma 
at finite temperatures and densities in the investigation 
of chiral symmetry restoration. With increasing tem- 
perature and density sigma changes its character from a 
resonance with large mass to a bound state with small 
mass. This medium dependence of sigma mass makes im- 
portant sense not only in thermodynamics [8], but also 
in particle yields. Around the critical point of the chi- 
ral transition sigma is one of the most numerous species, 
since it is nearly massless. If the transition happens at 
the freeze-out these light sigmas will lead to large event- 
by-event fluctuations [9] of final state pions due to the 
cr-exchange in the pion thermodynamic potential [10]. 

A direct signature of these light sigmas in hot and 
dense matter is a nonthermal excess [10] of pions due 
to the decay process a — > 2tt. With the expansion of 
the hadronic system produced in relativistic heavy ion 
collisions, the in-medium sigma mass rises towards its 
vacuum value and eventually exceeds the 7T7t threshold. 
As the <77T7r coupling is large, the decay proceeds rapidly. 
When this process occurs before freeze-out, the produced 
pions will be thermalized in the heat bath. However, 
when the decay happens after freeze-out, the generated 



pions do not have a chance to thermalize. Thus, it is ex- 
pected [10,11] that the resulting pion spectrum will have 
a nonthermal enhancement at low momentum. 

The dependence of non-strangle hadron masses on tem- 
perature T and baryon density n can be well parameter- 
ized as [12] 

m(T, n) ■ ' t \ 1 




where T c is the critical temperature of chiral phase transi- 
tion, and uq the baryon density of normal nuclear matter. 
From this parameterization, the temperature and density 
effects are quite different: The temperature dependence 
is very weak at low temperatures (T/T c < 1/2), while the 
hadron masses decrease linearly in density in the whole 
density region. This difference between temperature and 
density effects is reflected in the pion spectra in high en- 
ergy nuclear collisions. For ultra-relativistic heavy ion 
collisions where the central region is almost baryon free 
and the temperature effect is dominant, the sigma and 
pion masses are nearly constants when the temperature 
of the system is much below the critical one. There- 
fore, the pions produced from nonthermal sigma decay 
after fr eeze-out will have almost the same momentum 

Po = \J m "^'° - — rri*(Q, 0) in the vacuum. Aside from an 
overall enhancement of pions due to the thermal sigma 
decay, there will be also a peak in the pion spectra due 
to the nonthermal sigma decay. This peak will disappear 
partly when the density effect can not be neglected, since 
in this case the sigma mass varies significantly even after 
the freeze-out. Therefore, the sharp peak in the spec- 
trum will be replaced by a smooth low-momentum pion 
enhancement. 

We now illustrate the above statement by using a sim- 
ple statistic model. The measured final-state pions in- 
clude the direct pions emitted from the fireball and the 
pions from sigma decay, 



N = N dir + N de 



(2) 



Ndir is determined at freeze-out, and N de c is twice the 
number of sigmas at chiral phase transition where sigma 
has minimum mass and is then most numerous, 



Ndir = V f n^(T f ) , 

N dec = 2N a (T c ) = 2V c n a (T c ) 



(3) 



where Tf and Vf are the temperature and volume of the 
system at freeze-out, T c and V c the temperature and vol- 
ume at chiral phase transition, n w (T/) and n CT (T c ) the 
pion and sigma number densities at Tf and T c , 
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energies 



e v (T f ) = y/ml(T f )+p2 and e CT (T c ) = ^m2(r c )+p2. 
In central region of relativistic heavy-ion collisions where 
there is almost no net baryon and one considers the tem- 
perature effect only, the Bjorken's scaling hydrodynamics 
[13] is often used to describe the space-time evolution of 
the produced hot system. In this scenario the tempera- 
ture Tf , volume Vf and proper time r/ at freeze-out and 
the corresponding quantities T c , V c and t c at chiral phase 
transition have simple relations 



T f = T r [ — 
V f = V c T -L. 



(5) 



For central A — A collisions with impact parameter 
6 = 0, the number of participant nucleons can be taken 
as (N par t) = 2 A and the volume at freeze-out can be 
assumed to be proportional to A, Vf — Avf. The pseu- 
dorapidity density of charged pions per participant pair 
in the mid-rapidity range \r]\ < 1 is then written as, 
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where n„(Tf) is the pion number density in the angle 
interval 2arctane _1 < 6 < 2arctane, 
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and the factor of 2/3 is the charged to total pion ratio. 

The pseudorapidity density (6) depends on the vol- 
ume parameter Vf and the two temperature scales Tf 
and Td which characterize the collective effect of the sys- 
tem at freeze-out and at chiral transition. The volume 
parameter Vf is centrally colliding energy dependent, it 
can be fixed by experimental data. From the lattice 
simulations of QCD which make measurements of the 
Polyakov loop and the scalar quark density, the decon- 
fincment and chiral transitions coincide at a temperature 
of about T c = 170 MeV [14]. In the QCD phase diagram 
[15] the thermal freeze-out temperature in high energy 
nuclear collisions is about Tf = 120 MeV. As is well 
known, the pion mass is approximately temperature in- 
dependent in the whole chiral breaking phase, we can 
take (Tf) = m 7r (0) = 140 MeV at freeze-out and 
m a (T c ) = 2m 7r (T c ) = 2m n (0) at chiral phase transition. 



With the given Tf and T c , the ratio of the pions produced 
by sigma decay to the direct pions is 



n a (T c 



52% 



(8) 
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this shows a significant influence of the chiral phase tran- 
sition on the final state pions in relativistic heavy-ion 
collisions. 

In above study of pion multiplicity we did not distin- 
guish the nonthermal sigma decay from thermal sigma 
decay. As we discussed qualitatively, the thermal and 
nonthermal pions have very different momentum spec- 
tra. The pions produced by sigma decay in the tempera- 
ture region Tf < T < T c have time to thermalize before 
freeze-out, they lead to an enhancement of thermal pions, 
while in the interval T < Tf the produced pions do not 
get a chance to thermalize, they have almost the same 

momentum p = \J m "^'°^ - — m^(0, 0) in the rest frame 
of sigma due to the weak temperature dependence of non- 
strangle hadron masses at low temperatures. Therefore, 
if we take a narrow momentum window around po, only 
a small part of the thermal pions including the direct 
pions and thermal pions from sigma decay distribute in 
the window, but all the nonthermal pions are deposited 
it. The contribution from sigma decay to final state pions 
is therefore greatly amplified in the window. This strong 
pion enhancement can be taken as a signature of chiral 
phase transition in relativistic heavy- ion collisions. 

Considering a symmetric window Ap around po and 
taking the thermal fraction of the pions produced by 
sigma decay to be (3, the pseudorapidity density in the 
mid-rapidity region and in the momentum window is 
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with the pion number density 



^(T f ) 



d 3 p 



(9) 



(10) 



AO.Ap JaBAp ( 2?r ) 3 et " /Tf - 1 

in the window. 

From (9) the pion enhancement factor due to thermal 
and nontheral sigma decay in the window is 
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it depends on the thermal fraction (3 and the momentum 
window width Ap. 

In our numerical calculation we take the sigma mass in 
the vacuum m CT (0, 0) = 600 MeV, the corresponding mo- 
mentum of the nonthermal pions is then po = 265 MeV. 
For a symmetric momentum window around p the en- 
hancement factor rAp is shown in Fig.(l) as a function 
of the window width for three values of the thermal frac- 
tion (3 and in Fig. (2) as a function of the thermal fraction 
(3 for three values of the window width r& p . The en- 
hancement drops down rapidly with increasing window 
width and becomes saturated when the width is larger 
than about 0.4 GeV. It is inversely proportional to the 
thermal fraction. In the limit of wide window or large 
thermal fraction, the nonthermal sigma decay can be ne- 
glected and the enhancement in the window approaches 
the total enhancement (8). It is clear to see that the pion 
number in a small momentum window increases dramat- 
ically due to nonthermal sigma decay. Even for (3 = 0.9, 
namely only 10% of the sigmas decay after the freeze-out 
of the system, the enhancement in a narrow window is 
still significant. 



cial en hancement in a small momentum window around 
Po = \J ■^ 2 - — m% ~ 260 MeV. This anomalous enhance- 
ment comes from the temperature effect which does not 
affect the sigma mass at low temperatures. For density 
effect, the sigma mass drops down rapidly even at low 
densities, the pions from sigma decay after the freeze-out 
distribute also in a wide momentum region, then the pion 
enhancement becomes weak in the window. Therefore, 
in central region of relativistic heavy-ion collisions where 
there is almost no density effect, the pion enhancement 
in a small window around the nonthermal momentum p 
may be considered as a signature of the chiral phase tran- 
sition. Although the sigma width and the sigma momen- 
tum which are not taken into account in our estimation 
will smooth the sharp peak around po, we expect that 
the dynamics of chiral phase transition (here via temper- 
ature dependent masses) plays still an essential role in 
the pion spectra. 

acknowledgments: This work was supported in part 
by the NSFC under contract numbers 19925519 and 
10135030, and by the Major State Basic Research Devel- 
opment Program under contract number G2000077407. 




o l . . ■ ■ 1 

0.2 0.4 0.6 0.8 1 

Ap (GeV) 

FIG. 1. The window width dependence of the enhancement 
factor (11) for f3 = 0.1,0.5 and 0.9. 
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FIG. 2. The thermal fraction dependence of the enhance- 
ment factor (11) for Ap = 0.1,0.3 and 0.5 GeV. 

In summary, we have investigated the effect of chiral 
phase transition on final state pion spectra. The sigma 
enhancement at chiral phase transition leads to a signif- 
icant enhancement of pions due to the sigma decay in 
chiral symmetry breaking phase. Especially, the sigma 
decay after the freeze-out of the system results in a cru- 
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